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ABSTRACT 

In the early phase of the Solar System evolution, while the outer planets migrated due 
to their interaction with a planetesimal disk, Ju piter may have crossed the 2:1 mean 
motion resonance with Saturn. It is well known (jMorbidelli et al.ll2005f ) that this dy- 
namical event has profound consequences on the evolution of an alleged initial Trojan 
population of Jupiter. In this paper, we analyse in details the dynamics of Jupiter 
Trojans during the resonance crossing. We find that orbital instability is not confined 
to the central 2:1 resonance region but occurs in a more extended region where a sec- 
ular and secondary resonances perturb the Trojan orbits while the planets approach, 
cross and leave the 2:1 resonance. In addition, Jupiter and Saturn are locked after the 
resonance crossing in an apsidal corotation which has an additional destabilizing effect 
on Trojans. The synergy of the secular resonance, secondary resonances and apsidal 
corotation is needed to fully remove an initial Trojan population. New Trojans can be 
temporarily captured from the planetesimal disk while Jupiter crosses this extended 
instability region. After the disappearence of major secondary resonances, the secular 
resonance and the break of the apsidal corotation, the temporarily captured Trojans 
are locked and can remain stable over long timescales. 

Key words: celestial mechanics - asteroids. 



1 INTRODUCTION 

According to a widely accepted scenario, the outer planets 
of the solar system are embedded in a gas-free disk of plan- 
etesimals in the last stage of planetary formation. Gravita- 
tional interactions between planets and planetesimals dom- 
inate the dynamical evolution of the disk. Planetesimals are 
scattered by the planets in a chaotic manner. Orbital angular 
momentum and energy are exchanged resulting in planetary 
migration eventually with significant orbital changes. This 
phenomenon, described extensively for the first time by Fer- 
nandez and Ip (1984), was invoked then to suggest Pluto's 
capture in a 3/2 resonanc e by Neptune during its outward 
migration l|Malhotralll993h . 

At what heliocentric distances did the planets form and 
how far did they migrate? Two major models are proposed: 
in the first model investigated by Fernandez and Ip (1984) 
and applied by Malhotra (1993), initial planetary orbits are 
widely spaced between 5.2 and a bout 25 AU. A rec ently 
proposed model, the NICE model (|Tsiganis et al, 2005), as- 
sumes an initially closely spaced distribution between 5.3 
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and 17 AU. In the latter model, Uranus and Neptune ex- 
change their orbits during migration. 

The two migration models differ mainly in the assump- 
tion of the initial semimajor axes of the planets. The driv- 
ing mechanism for planetary migration is the same: planets 
scatter planetesimals in-and outwards. Inwards scattering 
moves a planet outwards while it moves inwards when a 
planet scatters a planetesimal outwards. In a closed system 
in equilibrium without loss of planetesimals where planets 
scatter planetesimals in- and outwards, no significant mi- 
gration would occur. In an open system with loss of plan- 
etesimals, on the other hand, significant planetary migration 
is possible. Jupiter plays a crucial role since it ejects easily 
planetesimals received from the other planets out of the so- 
lar system. As a consequence, Jupiter migrates towards the 
Sun while the other three planets migrate outwards. Migra- 
tion is halted when the outermost planet reaches the edge 
of the planetesimal disk and when most of the planetesimals 
scattered between the planets are removed. 

Jupiter and Saturn cross in the NICE model the 2:1 
mean motion resonance (MMR) soon after planetary mi- 
gration has started. As a consequence, the eccentricities of 
both planetary orbits increase. Saturn approaches the orbit 
of the third planet which is excited and which, therefore, 
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has close approaches with the fourth planet. The third and 
fourth planet may exchange orbits which moves the third 
planet rapidly towards 20 AU deep inside the planetesimal 
disk surrounding in the beginning the four planets. Dynam- 
ical friction with planetesimals damps rapidly enough plan- 
etary eccentricities to avoid close encounters between the 
third and fourth planet which would result eventually in a 
destabilization of the outer planetary system. The planetary 
orbits separate due to migration and their eccentricities are 
damped to present values due to dynamical friction. 

In a scenario where Jupiter and Saturn cross the 2:1 
MMR, Jupiter Trojans are destabilized. The destabilization 
was first attribu ted to the particular pe rturbations of the 
2:1 MMR solely l|Michtchenko et al.ll200lh . Later, Morbidelli 
et al. (2005) attributed the destabilization to the 3/1 sec- 
ondary resonance between harmonics of the libration fre- 
quency of Trojan orbits and a critical argument of the 2:1 
MMR. This secondary resonance is very effective to remove 
Trojans in case of a very low migration speed. Within a 
frozen model without migration, all Trojans are removed 
l|Morbidelli et al.ll2005h on a timescale of about 1 Myr. One 
cannot exclude, however, that a considerable fraction of Tro- 
jans survives since each secondary resonance is quite narrow 
and Trojans may pass through. In this paper, we will show 
that due to the presence of a major secular resonance on 
both sides of the 2:1 MMR, original Trojans are removed due 
to the synergy between secondary and secular resonances in- 
dependently of the planet migration rate. In addition, the 
lock of Jupiter's and Saturn's orbits into apsidal corotation 
after the 2:1 MMR crossing significantly contributes to the 
destabilization until the locking is broken. 

While primordial Trojans are destabilized before, dur- 
ing and after the crossing of the 2:1 MMR, nearby planetes- 
imals can be temporarily trapped on Trojan orbits via the 
reverse chaotic path. As soon as Jupiter leaves the extended 
instability region, the latest captured Trojans remain locked 
on tadpole orbits for long timescales comparable to the age 
of the planetary system. Morbidelli et al. (2005) have shown 
that the orbital distribution of the observed Trojans corre- 
sponds to the orbital distribution of the captured Trojans. 

Temporary trap ping in co o rbital motion appears to oc- 
cur still at present. lEverhartl (|l973l) desc r ibed temporary 
captures in horseshoe orbits and Karlssor] (|2004 ) identified 
about 20 transitional objects in a sample of about 1200 Tro- 
jans. Candidates are Centaurs that can be trapped as Tro- 
jans for short periods of 10 4 — 10 5 yrs dHorner and Evansl 
120061 ). This shows that the stabl e region for Jupiter T rojans 
is surrounded by a chaotic layer (|Marzari et al.ll2003h where 
a population of temporary Trojans resides. At present, the 
stable and unstable regions are well separated and an ob- 
ject residing in the transient population cannot become 
a permanent Jupiter Trojan without the help of a non- 
conservative process. There are some slow diffusion gates 
fro m the stable to the unstabl e region like those identified 
by l|Robutel and Gabern|[2006h related to commensurabili- 
ties between the secular frequency of the Trojan perihelion 
longitude and the frequency of the Great Inequality (2:5 al- 
most resonance between the present Jupiter and Saturn). 
However, it is very unlikely that a transient Trojan can fol- 
low in the reverse sense these paths to became a permanent 
Trojan. 

We describe in section 2 the major perturbations acting 



near the 2:1 MMR on Jupiter Trojans in the early phase of 
the NICE migration model. Section 3 is devoted to the syn- 
ergy of these major perturbations leading to a total loss of a 
possible initial Trojan population independent of migration 
rates. In section 4 we show that perturbations in the central 
Jupiter-Saturn 2:1 MMR region, where at least one of the 
resonant arguments librates, do not lead to global instability 
even in a frozen model as suggested by Michtchenko et al. 
(2001). A Frequency Map Analysis reveals extended stable 
regions. 



2 SOURCES OF INSTABILITY FOR JUPITER 
TROJANS BEFORE, DURING, AND AFTER 
THE 2:1 MMR CROSSING 

In this Section we describe the sweeping of resonances 
through the Trojan region before and after the 2:1 MMR 
crossing of Jupiter and Saturn. 



2.1 The numerical models 

The goal of our numerical modellig is to explore the stability 
of Trojan orbits during the migration of Jupiter and Saturn 
through the 2:1 MMR. The migration rates of the two plan- 
ets have to be computed within a model that includes all 
the outer pl anets and a disk of massive planetesimals as 
described in l|Tsiganis et al.| [2005). For this reason we have 
first reproduced the dyna mical evolution of the outer planets 
using the same model of (jTsiganis et al.ll2005l) and adopt- 
ing t he same SYMBA5 numerical algorithm (|Duncan et al.l 
1 1998T ) . It is a symplectic integrator that models the gravi- 
tational interactions among planets, the gravitational forces 
exerted by the planets on planetesimals and vice versa. The 
gravitational interactions among planetesimals is omitted in 
order to gain computing time. SYMBA5 is particularly de- 
signed to handle close encounters among planetesimals and 
planets, the main mechanism driving the migration of the 
outer planets. Using the starting co nditions for the plan- 
ets described in (|Tsiganis et al.|[2005l ) for relatively slow mi- 
gration, we performed a numerical simulat ion that matches 
closely that shown in (jTsiganis et al.|[2005l ). Hereinafter, we 
refer to this simulation as PPS (Planets and Planetesimals 
Simulation). The four outer planets are started on closely 
packed, almost circular and coplanar orbits. The semimajor 
axes a of Jupiter and Saturn are 5.45 and 8.50, respectively, 
so that they will cross the 2:1 MMR during their migra- 
tion. Orbital eccentricitie s e and inclinations i are equal to 
0.001 at start. Following (|Tsiganis et alj|2005l ), we use 4500 
massive planetesimals to produce the migration of the four 
planets. 

In Fig. [1] we show the semimajor axis and eccentricity 
of Jupiter as obtained in our PPS simulation. Before the 2:1 
MMR crossing, Jupiter's eccentricity is equal on average to 
0.01 in spite of its small starting value. This is due to the 
forced component of Saturn which grows while approaching 
the resonance location. The 2:1 MMR crossing is marked 
by a sudden jump in eccentricity related to the separatrix 
crossing and by large oscillations in semimajor axis. After 
the crossing, the eccentricity is slowly damped down while 
the planet continues to migrate towards its present location. 

From the time series of the orbital elements of both 
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Figure 1. Orbital evolution of Jupiter semimajor axis and eccen- 
tricity during the migration of the planets driven by planetesimal 
scatterin g (PPS simulat i on). T he outcome is very similar to that 
shown in Tsiganis et al. (20051). We focus on the resonance cross- 
ing that occurs after about 1.8 Myr from the beginning of the 
simulation. 



Jupiter and Saturn, computed within the PPS simulation, 
we can derive the planet migration rate da/dt of the semi- 
major axis and the eccentricity damping rate de/dt and pro- 
duce a synthetic model. In this model the effect of the plan- 
etesimal scattering is simulated by adding analytically the 
da/dt and de/dt terms to the equations of motion of the 
planets. Such an approach was exploited to model the ef- 
fect of circum s tellar disks on exopla n ets, f or instance by 
iLee and Peak] (|2002l ) and iKlev etaH (|2004l ). The authors 
used analytic expressions to estimate the changes in a and 
e due to the interactions with the disk when advancing the 
planets from time tj to ti+\. We follo w the formalism out - 
lined in the appendix of the paper bv lLee and Peald {2002) 
and, to model the migration of planets, we introduce a da/dt 
and de/dt in the SYMBA5 integrator to produce the migra- 
tion and neglect all the massive planetesimals. We concen- 
trate on the orbital evolution of Jupiter and Saturn since 
thay are responsible for the stability or instability of Jupiter 
Trojans. Uranus and Neptune are needed in the PPS model 
in order to transport the planetesimals responsible for the 
migration of the outer planetary system. However, by a se- 
ries of numerical tests, we have verified that their influence 
on the Trojan orbits of Jupiter is negligible compared to 
that of Jupiter itself and Saturn. 

In the synthetic model we must account for the fact 
that the migration of Jupiter and Saturn caused by plan- 



etesimal encounters is linear only over a limited amount of 
time and not over the whole migration period. The number 
of planetesimals in planet crossing orbits is in fact declining 
causing a slow decrease of da/dt and de/dt. We, therefore, 
tune the synthetic integrator to the PPS run by using values 
of da/dt and de/dt that are derived from PPS at different 
times during the evolution of Jupiter and Saturn. In this 
way, the synthetic model accurately reproduces the evolu- 
tion of the planets during the 2:1 MMR and even after. 
Moreover, it retains all the dynamical features needed to 
analyse the stability of Jupiter Trojans. The initial values of 
da/dt and de/dt for Jupiter are —7.39 x 10 -9 AU/year and 
— 3.76x 10 -10 1/year, respectively. After 10 Myr these values 
have decreased to -4.05 x 10" 9 AU/year and —2.05 x 10~ 10 
1/year. For Saturn, the da/dt ranges from 2.23 x 10~ 8 to 
1.24 x 10~ 8 AU/year in 10 Myr while the de/dt goes from 
-2.80 x 10" 9 to -1.65 x 10~ 9 1/year. 

The main advantage of using the synthetic integrator is 
its speed. We can compute the orbital evolution of Jupiter 
and Saturn and of massless Trojans on a timescale at least 
100 times shorter than that required by a full model that 
includes the massive planetesimals (PPS-type model). CPU 
time is a critical issue since we have to explore the stability 
of Trojans in the phase space for different intervals of time 
during migration and in different dynamical configurations. 
In addition, with the synthetic model there is the possibility 
of easily changing the values of da/dt and de/dt which are 
strongly model dependent. We tested also a synthetic model 
based on the RADAU integrator and the results were in 
agreement with the SYMBA5 synthetic model. 

To identify possible resonances between the motion of 
the planets and that of Trojans we have to evaluate the ma- 
jor orbital frequencies of these bodies. However, the dynam- 
ical system evolves because of planetary migration and the 
frequencies change with time. To compute the value of these 
frequencies at a given instant of time we use frozen models. 
We extract the osculating orbital elements of the planets 
and Trojans at the required time and start a numerical in- 
tegration of the trajectories with the migration switched off 
(both da/dt and de/dt are set to 0). In this way, we com- 
pute a time series of orbital elements for the non-migrating 
planets and Trojans long enough to derive precise values of 
the frequencies. 

To compute initial orbital elements for Trojans at dif- 
ferent times during the evolution of the planetary system 
we select random initial conditions within a ring surround- 
ing the orbit of Jupiter. The semimajor axis of any puta- 
tive Trojan is selected in between 0.9 aj and 1.1 a.j where 
a.j is the semimajor axis of Jupiter. The eccentricity can 
be as large as 0.5 and the inclination extends up to 50°. 
The other orbital angles are selected at random between 
and 360°. Each set of initial conditions is integrated for 10 4 
years and if the critical argument A — Aj librates in this 
time period, a body with that set of initial conditions is in- 
cluded in the sample of virtual Jupiter Trojans. The choice 
of wide ranges in eccentricity and libration amplitude, some- 
what wider than the present ones, is dictated by the chaotic 
evolution of the orbital elements before, during and after 
the 2:1 MMR crossing. This chaotic evolution can drive a 
given Trojan orbit from a high eccentric orbit into an al- 
most circular one, and it can strongly reduce its libration 
amplitude. We cannot neglect at this stage orbits that are 
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Figure 2. Period corresponding to frequency g\ as a function of 
Jupiter's semimajor axis. 

unstable on the long term since they might be turned into 
stable ones during the dynamical evolution caused by the 
planetary migration. A body is considered to be ejected out 
of the swarm during its evolution when its critical argument 
no longer librates. 

2.2 Secular resonance with Jupiter 

The secular evolution of eccentricities and perihelion lon- 
gitudes of the Jupiter-Saturn system, as described by the 
Lagrange-Laplace averaged theory, is characterized by two 
major frequencies that we call g\ and t/2 following Murray 
and Dermott (1999). These frequencies are not constant dur- 
ing planetary migration since their values depend on the 
semimajor axes of the two planets through the Laplace co- 
efficients. The linear Lagrange-Laplace theory has an ana- 
lytical solution that allows to compute both gi and gi as a 
function of planetary orbital elements. However, this solu- 
tion fails in proximity of the 2:1 MMR and we resort to a 
full numerical approach to compute the two frequencies dur- 
ing planetary migration. We 'freeze' the dynamical system 
at different stages of migration (frozen model) and we esti- 
mate both g\ and g2 from the time series of the non-singular 
variables h and k of Jupiter over 1 x 10 6 yrs. As usual, we de- 
fine these variables by h = e*cos(tu) and k — e*sin(tu). For 
the computation of precise values for the two frequencies we 
use t he so-called Modif i ed Fourier Tra nsform (MFT) anal- 
ysis (|Laskar et al.l[l9 92; Laskar 1993a, b), which we had al- 
ready applied to study the s tability properties of the present 
Jupiter Trojan population dMarzari et al.l [20031 ). gi and <?2 
are by far the frequencies with the largest amplitude com- 
puted from the MFT. 

One of the two frequencies sweeps through the Trojan 
region during the migration of the planets reaching values 
typical of the proper frequency g of Jupiter Trojans. We call 
this frequency g\ while the other frequency, <?2, has a longer 
period and does not influence the Trojan motion. When g\ 
is equal or very close to g a secular resonance is established. 
Fig. [5] and Fig. [3] show the behaviour of g\ as a function of 
Jupiter's semimajor axis and of time during migration, re- 
spectively. The period corresponding to g\ decreases while 
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Figure 3. Period corresponding to frequency g\ as a function of 
time. 

Jupiter and Saturn approach the 2:1 MMR and it rises back 
after the 2:1 MMR. Fig. [4] shows for comparison the ratio 
of the orbital periods between Saturn, Ps, and Jupiter, Pj. 
When Jupiter and Saturn approach, cross and leave the 2:1 
MMR, <?i sweeps through the Trojan phase space causing 
strong perturbations that lead mostly to instability. Libra- 
tion amplitudes and/or eccentricites of Trojans are increased 
resulting in close encounters with Jupiter. Due to the func- 
tional dependence of a on th e proper elements of the Trojan 
orbits (|Marzari et al.ll2003l ) the secular resonance appears 
first at high inclinations, moves then down to low incli- 
nations when the planets reach the 2:1 MMR, and finally 
climbs back to high inclinations after resonance crossing. 
This behaviour will be described in more detail in Section 
3.1. Fig. [S] shows a power spectrum of the complex signal 
h + ik for a Trojan orbit. The frequency g\ approaches g 
when the planets migrate towards the 2:1 MMR leading to 
resonant perturbations. The frequency gi does not change 
much and remains far from g. The Trojan becomes unsta- 
ble just after the third instant of time as shown in Fig. [5] 
before the 2:1 MMR is crossed (t=1.5 Myr). When it falls 
inside the g = g\ resonance its orbit is in fact destabilized 
on a short timescale by a fast change in eccentricity and 
libration amplitude. 

However, the delicate dynamical equilibrium of the Tro- 
jan motion is perturbed even when g\ is only close to 
g, outside the secular resonance borders. The term pro- 
portional to g — gi in the disturbing function is dynam- 
ically important generating a chaotic evolution of Trojan 
orbits even if on a longer timescale compared to those cases 
falling into the resonance. A similar effect was observed for 
Uranus Trojans whose diffusion speed in the phase space is 
strongly increased, leading to chaotic motion, in proximity 
of the fundamental f requencies <js and 57 of the solar system 
l|Marzari et alj|2003h . When gi leaves the Trojan region af- 
ter the 2:1 MMR, it remains anyway close to g for a long 
time persisting as a source of instability. Moreover, after 
the 2:1 MMR, Jupiter and Saturn are locked in an apsidal 
resonance that enhances the strength of the g — gi term by 
coupling the perturbations of Jupiter to those of Saturn (see 
next section). 
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Figure 4. Ratio between Saturn's (Ps) and Jupiter's (Pj) orbital 
periods as a function of time. 
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Figure 5. Power spectrum of h and k for a Trojan while Jupiter 
and Saturn approach the 2:1 MMR. The secular frequencies g\ 
and <?2 of the Jupiter-Saturn system and the Trojan's proper fre- 
quency g are given at three different times. Frequency gi moves 
towards g destabilizing the Trojan's orbit at 1.5 Myr. 



A change in the initial values of Jupiter and Saturn 
in the migration model would move the location of the 2:1 
MMR and the corresponding values of the semimajor axes of 
both Jupiter and Saturn at the crossing. However, this does 
not alter the effect of the secular resonance on the stability of 
Trojans. The resonance sweeping occurs anyway since gi and 
<?2 depend on the semimajor axis of Jupiter a.j, according 
to the Lagrange-Laplace averaged theory, in the same way 
as the frequenc y g depend s on a.j following Erdi's theory of 
Trojan motion (|Erdill 19791 ). 



2.3 Secondary resonances with harmonics of the 
'2:1 Great Inequality' 

There are two independent critical resonance arguments for 
the 2:1 MMR of Jupiter and Saturn: f?i = Aj-2A s +otj and 
#2 = \ j — 2Xs+vos, where A and cc denote respectively mean 
longitude and longitude of perihelion. Either one of the two 
critical arguments librates while the other circulates or both 



o 



CD 



1000 
500 


-500 
-1000 



I I 


i i 






0^ 




,e— e 






3:1 2:1^ f 






I 


5 1 



0.0e+00 1.0e+06 2.0e+06 3.0e+06 
Migration time (yr) 

Figure 6. Evolution of the circulation period Tf of f?i as a func- 
tion of time. The resonance crossing is marked by a discontinuity 
in the period of 8\ . 



critical arguments librate simultaneously. In the latter case, 
the difference between the two critical arguments 02 — 9\ = 
Avj also librates. This means that Jupiter and Saturn are 
in apsidal corotation. 

While Jupiter migrates towards the Sun and Saturn in 
opposite direction, both 9\ and 62 circulate prograde be- 
fore the 2:1 MMR and retrograde after. The frequency of 0i 
and 6*2 may become commensurable with the libration fre- 
quency of the critical argument of Jupiter Trojans. This is 
the case of a secondary resonance which was investigated by 
Kortenkamp et al.(2004) for Neptune Trojans. The authors 
found that a Neptune Trojan in a secondary resonance can 
enhance significantly its libration amplitude possibly lead- 
ing in some cases to instability. The importance of secondary 
resonances for Jupiter Trojans in the frame of the NICE 
model was recognized by Morbidelli et al. (2005). Secondary 
resonances can be encountered before and after the cross- 
ing of the 2:1 MMR. In a frozen model without migration, 
the 3:1 secondary resonance after the 2:1 MMR removes all 
Trojans on a timescale of 1 Myr while the 2:1 secondary res- 
onance removes 70% of them. In a migration model, these 
removal rates can be significantly less if the secondary reso- 
nances are crossed rapidly. 

In Fig. [5] we show the period Tf of the frequency / 
of circulation of 81 as a function of the semimajor axis of 
Jupiter during migration. Different secondary resonances are 
crossed. Crossing, however, is fast, in particular after the 
2:1 MMR. In proximity of the 4:1 secondary resonance, for 
example, the period of 81 changes by approximately 20% 
in only 3 x 10 4 yrs. In Fig. [7] we illustrate with a shaded 
stripe the frequency interval of / (translated into periods) 
for which there is a 4:1 (lower shaded stripe), 3:1 (middle 
shaded stripe) and a 3:1 commensurability (upper shaded 
stripe) with the libration frequency of a Trojan swarm. We 
consider Trojans up to 50° in inclination and up to 0.35 
in eccentricity corresponding to libration periods roughly 
ranging from 145 and 190 yrs. The sweeping appears to 
be fast, in particular for the 4:1 and 3:1 secondary reso- 
nances, taking also into account that any individual Trojan 
will be affected only by a fraction of the time spent by / 
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Figure 8. Jumps in the libration amplitude of a Trojan while 
crossing the secondary resonances 2:1 near 1 Myr and 3:1 near 2 
Myr before the 2:1 MMR. The Trojan enters at the end near 2.5 
Myr the secular resonance g\ and is destabilized. 



to cover the entire shaded region. It is worthy to note here 
that the migrat ion speed is relatively low within the differen t 
NICE models dTsiganis et al.ll200a iMorbidelli et alj 120051 ). 
A faster migration would further reduce the relevance of sec- 
ondary resonances in the destabilization of Jupiter Trojans 
during the 2:1 MMR. 

As for the secular resonance, the crossing of the sec- 
ondary resonances occurs before and after the 2:1 MMR. 
However, there is a substantial difference between the two 
dynamical configurations. Before the 2:1 MMR the sec- 
ondary resonance sweeping causes sharp jumps in libration 
amplitude and eccentricity that in most cases do not fully 
destabilize the Trojan orbit. As shown in Fig. [S] the cross- 
ing of the 2:1 secondary resonance at t ~ lMyr reduces 
the libration amplitude increasing the stability of the orbit. 
When the 3:1 secondary resonance is encountered later, the 
initial libration amplitude is restored. The Trojan orbit be- 
comes finally unstable when it crosses the secular resonance 
with g\. Of course, for librators with large amplitudes, the 
perturbations of the secondary resonances may lead to a 
destabilization of the Trojan orbit. 

Totally different is the dynamical behaviour after the 
2:1 MMR. The secondary resonances are much more effec- 
tive in destabilizing Trojan orbits independent of their li- 
bration amplitude. The reason for the different efficiency of 
secondary resonances before and after the 2:1 MMR is due 
to two independent causes: 

• Immediately after the 2:1 MMR crossing, the eccen- 
tricity of Jupiter is on average higher. This reinforces pre- 
sumably secondary resonances and the secular resonance. 
We tested this hypothesis before the 2:1 MMR by numer- 
ically integrating the same Trojan orbits in a model with 
the eccentricity of Jupiter set to an average value 2.5 times 
higher compared to that of the reference model which ap- 
proximately corresponds to the average increase observed in 
simulations. Trojans are started in between the 2:1 and 3:1 
secondary resonance. Trojans surviving at least 2 x 10 4 years 
decreased by 33 % with respect to the Jovian low eccentric- 
ity case. On average, the lifetime was reduced by 22 % in 
the higher eccentric case. 



• After the 2:1 MMR crossing the planets are always 
locked in apsidal corotation, according to our simulations. 
This additional dynamical effect contributes to destabilize 
Trojan orbits. To estimate the effects of apsidal corotation 
we have used the same model described in the previous item 
(that with higher eccentricity) and forced apsidal corotation 
of the planets before the 2:1 MMR crossing by a convenient 
choice of the orbital angles of the planets. A comparison 
between apsidal and non-apsidal corotation after resonance 
crossing is not possible since the system finds always rapidly 
the apsidal corotation state. In the apsidal corotation model 
the number of surviving Trojans drops by about 23 %, com- 
pared to that without apsidal corotation, and the Trojan 
lifetime is shortened by 42 %. 

2.4 Effect of apsidal corotation between Jupiter 
and Saturn on the dynamics of Trojans 

After the 2:1 MMR crossing, Jupiter and Saturn are locked 
in apsidal corotation in all our simulations. In most cases, 
apses are antialigned with Avj — zuj — zus librating about 
180°. This apsidal corotation is broken much later. The pres- 
ence of apsidal corotation, as stated in the previous section, 
has significant consequences for the instability of Trojans: 

• It enhances the effects of the g—gi secular term since the 
frequencies of the precessional rates for the perihelia longi- 
tudes of Jupiter and Saturn are commensurable. In Fig.|5]we 
show the power spectrum of a Trojan started in between the 
2:1 MMR and the 3:1 secondary resonance. The two peaks 
corresponding to the gi and gi frequencies of the Jupiter 
and Saturn system are clearly visible and g\ is close to the 
proper frequency g. The peaks are much higher as compared 
to the power spectrum in Fig. [S]which is obtained before the 
2:1 MMR where Attj circulates. 

• The secondary resonances after the 2:1 MMR become 
very effective in destabilizing Trojan orbits due to the in- 
creased eccentricity of Jupiter as pointed out above. The 
coupling between the apsidal corotation and secondary res- 
onances causes a fast growth of the eccentricity and a cor- 
responding shift in the libration center of the Trojan tad- 
pole motion away from the Lagrangian points L4 and L5 
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Figure 9. Power spectrum of a Trojan trajectory after the 2:1 
MMR crossing. The apsidal corotation increases the strength of 
the secular term gi that leads to chaotic evolution. 




20000 



Figure 10. Orbital evolution of a Trojan after the 2:1 MMR 
crossing perturbed by the 4:1 secondary resonance. The libration 
center is shifted from 300° to a lower value due to the apsidal 
corotation between Jupiter and Saturn. 
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Figure 11. Evolution of two fictitious Trojan populations while 
Jupiter and Saturn approach the 2:1 MMR. Population 1 is 
started 2.3 Myr before the resonance crossing. It is eroded mainly 
by the secular resonance g\ . Secondary resonances affect solely li- 
brators with large amplitudes. Population 2 is started 6 X 10 5 
yrs before the 2:1 MMR. Somewhat more bodies survive since 
the sweeping gi resonance affects less high inclined Trojans. Be- 
hind the 2:1 MMR, these surviving bodies are removed by effects 
discussed below. The continuous line marks the 2:1 MMR. The 
dotted horizontal lines show the location of the secondary reso- 
nances over all the range of libration amplitude of the bodies in 
the two populations. The timespan covered by the secondary res- 
onances is shrinking because the migration is faster in proximity 
of the 2:1 MMR. 



3 SYNERGY BETWEEN SECONDARY 
RESONANCES, THE SECULAR 
RESONANCE Gi AND APSIDAL 
COROTATION 

In order to investigate the combined effects of the three iden- 
tified major perturbations, we start Trojan populations at 
different migration stages before and after the 2:1 MMR 
crossing. Simulations begin in between the major secondary 
resonances 4:1, 3:1, 2:1, 1:1, and right before and after the 
onset of apsidal corotation. 

The starting values are produced by generating ran- 
domly the initial orbital elements and checking for the crit- 
ical libration angle. Maximal starting inclinations and ec- 
centricities are taken somewhat larger than in the presently 
observed Trojan population of Jupiter because, as already 
anticipated in Section 2.1, the chaotic evolution in the prox- 
imity of the resonance crossings may reduce their values and 
lead to a stable tadpole orbit. 



<|Namouni et al.lll999T ). In Fig. QTj] we show an example for 
this shift and eccentricity increase in the 4:1 secondary res- 
onance after the 2:1 MMR crossing. The Trojan is destabi- 
lized at the end by a close encounter with Jupiter. Note that 
this behaviour is never observed before the 2:1 MMR cross- 
ing where there is no apsidal corotation and the eccentricity 
of Jupiter is lower. 



3.1 Fate of Trojans before the 2:1 MMR crossing 

When Jupiter and Saturn approach the 2:1 MMR, Trojans 
cross secondary resonances and, in particular, the g\ secular 
resonance. Secondary resonances before the 2:1 MMR are a 
very weak instability source and destabilize solely tadpole 
orbits with large libration amplitude. The secular resonance 
g%, on the other hand, may remove Trojans with any li- 
bration amplitude when it sweeps through the region. Even 
when a body is not exactly within the gi resonance but 
closeby it feels the perturbations of the g — gi term and it 
may be destabilized, even if on a longer timescale. Fig. [11] 
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illustrates the erosion of two initial Trojan populations start- 
ing at different times, 2.3 and respectively 0.6 Myr before 
the 2:1 MMR. The two populations are generated with the 
same random process described in Section 2.1 and, as a con- 
sequence, they are dynamically similar. It appears that the 
secular resonance and the secondary resonances are more ef- 
fective close to the 2:1 MMR where they destabilize more 
than 90% of the Trojan population. The location of the sec- 
ondary resonances are shown during the evolution of the 
populations as horizontal lines in Fig. 1111 Before the 2:1 
MMR resonance crossing, the sweeping of these resonances 
is slower than after the the 2:1 MMR when the migration 
of the planets is much faster. However, as discussed above, 
before the 2:1 MMR secondary resonances are significantly 
weaker because of the reduced eccentricity of Jupiter. We 
also recall that each individual Trojan has its own libra- 
tion period and it is affected by the secondary resonances 
only during a fraction of the time covered by the resonance 
sweeping. 

The secular resonance sweeping through the Trojan re- 
gion is illustrated in Fig. [12] and Fig. 1131 Fig. [12] shows the 
escape time as a function of initial inclination for popula- 
tion 2. The secular resonance at the beginning of the simula- 
tion destabilizes a large number of bodies with inclinations 
around 25°. They all leave the Trojan region in less than 
1 x 10 5 yr. The critical value of 25° for inclination is de- 
termined by the choice of the initial orbits of Jupiter and 
Saturn before the 2:1 MMR which in turn determines the 
value of g%. A smaller initial distance between the planets, 
like in case of population 1, would have destabilized most 
Trojans at a higher inclination. As the planets move towards 
the 2:1 MMR during their migration, the frequency g\ in- 
creases and perturbs Trojan orbits at a progressively lower 
inclination. According to Marzari et al.(2003), the proper 
frequency g is higher for low inclined Trojans. In our sample 
of Trojans, there are naturally unstable orbits since we do 
not exclude those trajectories with large libration amplitude 
and high eccentricity. Some may be injected deeper in the 
stable region after crossing secondary resonances. Most of 
them, however, escape on a short timescale and populate 
the figure at the lower edge of the y-axis. 

The high inclined Trojans that survive the 2:1 MMR 
crossing are destabilized when the frequency gi decreases 
again (its period grows) as shown in Figs. [2] and Figs. [3] 
If some high inclination Trojans survive somehow the first 
sweep of the secular resonance, either by chance or because 
the planets start their migration close to the resonance lo- 
cation as in Fig. 1121 they probably will be destabilized by 
the second resonance sweeping when the planets move away 
from the resonance. 

In Fig. [13] we illustrate the distribution in inclination 
and semimajor axis of the same Trojan swarm integrated 
in frozen models with the planets progressively approaching 
the 2:1 MMR. The empty stripe corresponds to the secular 
resonance destabilizing the orbits on a timescale of 1 x 10 5 yr 
in a frozen model. The sweeping proceeds towards lower in- 
clinations while the planets approach the 2:1 MMR, whereas 
it rises back after the crossing in a symmetric way. We like to 
emphasize that bodies in the secular resonance are destabi- 
lized on a short timescale. Trojans whose frequency is close 
to gi but are not within the resonance borders are perturbed 
by the term g — gi and have a slower chaotic behaviour. This 
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Figure 12. Escape time for population 2 Trojans of Fig. 1111 vs. 
initial orbital inclination. The secular resonance sweeps through 
the Trojan region starting from an inclination of about 25° down 
to low inclined orbits. The continuous line marks the 2:1 MMR. 




Figure 13. Sweeping of the secular resonance when the planets 
approach the 2:1 MMR. The orbits of the same initial Trojan 
swarm are integrated within 4 different frozen models with the 
planets progressively closer to the 2:1 MMR. 



explains why also high inclined Trojans in Fig. [T^] are slowly 
eroded away. It also accounts for the fact that both pop- 
ulation 1 and 2 are fully destabilized notwithstanding that 
population 1 is started much farther away from the 2:1 MMR 
crossing. 



3.2 Evolution of a Trojan population after the 2:1 
MMR crossing: from the 4:1 to the 1:1 
secondary resonance. 

The closer Jupiter and Saturn start their migration to the 
2:1 MMR, the more initial Trojans may survive due to the 
dependence of the destabilizing secular resonance on orbital 
inclination. After the 2:1 MMR, the surviving Trojans would 
encounter a second time the secular resonance and the sec- 
ondary resonances. Both the secular and the secondary reso- 
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nances are reinforced by apsidal corotation and by the larger 
eccentricity of Jupiter and the destabilization rate is signif- 
icantly higher. In Fig. [14] we show the evolution of two pop- 
ulations of Trojans after the 2:1 MMR. The population in 
the upper diagram starts in between the 4:1 and 3:1 sec- 
ondary resonance. The Trojans of the second population in 
the lower diagram have libration frequencies which place 
them in between the 3:1 and 2:1 secondary resonances. To 
understand the features of the two diagrams, we have to keep 
in mind that the libration frequency depends on several or- 
bital elements. A Trojan population with about the same 
semimajor axis is spread over a large range of their other 
orbital parameters. As a consequence, the Trojans cross the 
same secondary resonance at different times. Moreover, the 
proper frequency g depends mostly on the inclination and 
the effects of the secular resonance appear, as already noted 
above in the discussion of Fig. 1121 at different inclinations 
during the sweeping. According to Fig. 1141 orbits with low 
inclination are rapidly destabilized while the escape time 
grows significantly for inclinations higher than 15°. At high 
inclinations the secular resonance arrives at a later time and 
it takes longer to destabilize Trojans. Some of the orbits 
of the second population survive the 3:1 crossing but are 
ejected before reaching the 2:1 secondary resonance. Only a 
few Trojans get beyond. At low inclinations where the in- 
stability is very fast, we observe a rapid pumping up of ec- 
centricity and a corresponding shift in the libration center. 
The power spectrum shows also the vicinity of the secular 
resonance. We conclude that, after the secondary resonance 
crossing, no initial Trojan can survive due to the synergy 
between secondary and secular resonances. Their effects are 
enhanced by the apsidal corotation of the two planets. The 
chaotic trapping of new Trojans appears to be difficult at low 
inclinations where the instability is very fast on timescales 
of the order of a few 10 4 yrs, while it might be more effi- 
cient at higher inclinations where the slow instability might 
allow the formation of a steady state transient population 
of unstable Trojans. 

At the 1:1 secondary resonance, significantly farther 
away from the 2:1 MMR, a sharp jump in libration am- 
plitude D and eccentricity occurs for Trojan orbits. An ex- 
ample is given in Fig. [K] where D changes during the 1:1 
crossing. This resonance is weaker as compared to the previ- 
ously encountered secondary resonances and it does not fully 
destabilize tadpole orbits but it induces chaotic variations 
of D. After the 1:1 crossing, when the libration frequency 
is away from the frequency of either 9i or 62, the libration 
amplitude still shows an irregular behaviour. By inspecting 
the power spectrum of the h and k variables of the Trojan 
orbit we find that the secular frequency gi is still relevant 
with a peak about half the size of the proper one. It is still 
a source of slow chaotic diffusion for the Trojan orbit. 

3.3 Far away from secondary resonances: still 
chaotic changes of orbital elements 

As noted before, when the system gets beyond the 1:1 sec- 
ondary resonance both the eccentricity and libration ampli- 
tude of tadpole orbits show a slow chaotic evolution which is 
enhanced when the planets cross mutual higher order mean 
motion resonances. Fig. \W\ shows the evolution of a Trojan 
trajectory with initially small values of D ~ 20°, e p ~ 0.03 
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Figure 14. Escape time vs. initial inclination for two Trojan 
populations started after the 2:1 MMR crossing. The first pop- 
ulation (top) has libration frequencies encompassed between the 
4:1 and 3:1 secondary resonance. The second population (bottom) 
is started in between the 3:1 and the 2:1 secondary resonance. 
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Figure 15. Evolution of the critical argument of a Trojan orbit 
crossing the 1:1 secondary resonance. The crossing is marked by 
the first sharp jump in libration amplitude. Subsequently, far from 
the secondary resonance, the orbit stays chaotic since the dynam- 
ics is still influenced by the the secular frequency g± powered up 
by the apsidal resonance. 
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Figure 16. The critical argument of a Trojan orbit started far 
away from secondary resonances. The orbit is still mildly chaotic 
because of the presence of the secular frequency gi . When Jupiter 
and Saturn cross the 4:9 mean motion resonance, the libration 
amplitude of the orbit increases until ejection out of the Trojan 
region. 

and i p ~ 19°. This orbit would lie deeply in the stable re- 
gion for the present configuration of the planets and it is far 
from any significant secondary resonance. The secular fre- 
quency gi appears to be still somehow relevant for the sta- 
bility of the Trojans causing moderate libration amplitude 
variations. However, the orbit is finally destabilized during 
the crossing of a 4:9 mean motion resonance between Jupiter 
and Saturn. The large libration amplitude increase begin- 
ning at ~ 5 x 10 6 yr leads to a fast destabilization of the 
tadpole orbit. Significant changes of the proper elements are 
still possible at this stage and the door for chaotic trapping 
is still open. 

After the apsidal corotation is broken, and the weak- 
ened secular frequency g\ has moved farther away from g, 
Trojan orbits are finally stable on a long timescale with no 
detectable variations of the libration amplitude. The door 
for chaotic capture is closed and the Trojan population ap- 
proaches its present configurations with no other significant 
remixing of proper elements. 



at least 10 5 yrs. Instability for this restricted starting region 
in phase space usually doe s not set on before 1 Myr. S imi- 
lar results were obtained b y lNesvornv and Donesl (|2002h and 
iMarzari and Scholll (|2002l , in static models were the planets 
were moving on fixed orbits. We will perform here a more 
detailed analysis of the stability of Trojans when Jupiter and 
Saturn are locked in the 2:1 MMR by using the FMA (Fre- 
quency Map Analysis) as described in Marzari et al.(2003). 
The semimajor axes of Jupiter and Saturn correspond to 
values of the NICE model. Migration is switched off, so that 
the planets do not leave the resonance (frozen model). As 
pointed out above, the two resonance variables 8i and 82 
may both librate (apsidal corotat ion), or only one m ay li- 
brate while the other circulates l|Marzari et al.l [20061 ). We 
applied the FMA analysis for both cases. Our results show 
that the stability of Trojans depends strongly on their initial 
conditions and on the behaviour of the two resonance argu- 
ments. The upper diagram in Fig. [T7] represents a diffusion 
portrait for Trojan orbits with Jupiter and Saturn in apsidal 
corotation. Corotation is possible around 0° or 180°. Since 
we obtain in most migration models corotation around 180° , 
we use this alignment mode for producing the diagram. Ex- 
tended stability regions appear between medium and high 
inclinations and for a large range of values for libration am- 
plitudes D of Trojans. Empty regions in the plot indicate 
instability times shorter than 1 Myr. The most stable region 
(the red one) has values for diffu sion speed comparab le to 
those of present Jupiter Trojans |Marzari et alj 2003) sug- 
gesting that bodies can survive for a long interval of time 
of the order of some Gyrs. For bodies with higher diffusion 
speed we still expect lifetimes of the order 10 7 — 10 8 yrs. The 
stable region extends down to low inclinations with libration 
amplitudes of about D ~ 60° where we found stability with 
different integrators, contrary to Michtchenko et al.(2001). 

In the lower diagram of Fig. [17] we consider a different 
dynamical state for the two planets in resonance. Only one 
of the two critical arguments librates. Consequently, A-cu 
circulates. The stability area is more extended in this case 
and orbits with low inclination can be found at low values 
of libration amplitude D. These results reinforce the idea 
that corotation contributes significantly to reduce dynamical 
lifetimes of Trojans. 



4 WHEN JUPITER AND SATURN ARE IN 
THE 2:1 MMR 

The instability of Jupiter Trojans with t he planets in the 2:1 
MMR with Saturn was investigated by iMichtchenko et all 
l|200lf ) in a frozen model. Jupiter and Saturn are in apsi- 
dal corotation. Trojan starting values are confined to incli- 
nations of 5°, zn — vj.j — 60°, A — A,/ = 60°, and eccen- 
tricity lower than 0.3. Using a RADAU integrator, the au- 
thors find instability over a very short time scale of about 
10 4 yrs. This indicates that if the migration of Jupiter and 
Saturn was very slow, a temporary capture of the plan- 
ets in the 2:1 MMR might have led to global instability 
of Trojans. However, when we performed numerical simu- 
lations of Trojan orbits in a frozen model like Michtchenko 
et al.(2001) we did not find short term instability. Using 
their semimajor axes for Jupiter and Saturn and confining 
Trojans to their starting region, and using also a RADAU 
integrator, we found a large number of stable Trojans over 



5 CONCLUSIONS 

We investigate the depletion of an alleged initial Jupiter Tro- 
jan population in the frame of the NICE model describing 
the early migration phase of the outer planets during which 
Jupiter and Saturn may have crossed their mutual 2:1 MMR. 
The loss of an initial population, possibly trapped during the 
growth of the planet, is due to the synergy of three different 
effects: 

• A secular resonance with the frequency g\, one of the 
two frequencies that, according to the Lagrange-Laplace 
theory, determine the the secular evolution of the eccen- 
tricity and perihelion longitude of Jupiter and Saturn. 

• Secondary resonances due to commensurabilities be- 
tween a critical resonance argument of the 2:1 MMR and 
the libration frequency of the critical argument of the Tro- 
jan orbits. 
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apsidal corotation of Jupiter and Saturn's orbit. When the 
frequency g\ moves out of the Trojan region but is still bor- 
dering it, the secular term g — gi is strong enough to perturb 
the Tojan motion causing instability on a relatively longer 
timescale. While Trojans are removed, new Trojans can be 
captured by the reverse chaotic path from the surrounding 
planetesimal population which drives planetary migration. 
The newly captured Trojans might be lost again until the 
secular resonance, secondary resonances and higher order 
mean motion resonances between Jupiter and Saturn disap- 
pear. 

The center of the 2:1 MMR, where at least one of the 
critical resonance arguments librates, is not particularly ef- 
fective in destabilizing Jupiter Trojans. Its effect is much 
weaker as compared to the secular resonance g\ and the sec- 
ondary resonances after the 2:1 MMR crossing. When the 
planets are steadily locked in resonance we find extended 
stability regions in the phase space of Trojan orbits. 
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Figure 17. Diffusion portraits of Trojan orbits for the 2:1 MMR. 
In the upper diagram both critical resonant arguments librate and 
the planets are in corotation. In the lower diagram, where a wider 
region has lower diffusion speeds, marked by red squares, only one 
critical argument librates. 

• Jupiter and Saturn's apsidal corotation after the 2:1 
MMR crossing 

While the planets approach the 2:1 MMR, the secular 
resonance gi sweeps through the Trojan region. It appears 
first at high inclinations and it moves down to almost zero 
degrees when the planets reach the center of the 2:1 MMR. 
It moves up again at higher inclinations after the resonance 
crossing, sweeping for a second time the Trojan region. Also 
secondary resonances appear before and after the 2:1 MMR 
crossing but they sweep across the Trojan region at a faster 
rate, in particular after the 2:1 MMR crossing. Before the 
2:1 MMR crossing, secondary resonances remove very few 
Trojans while they are more effective after the crossing be- 
cause of the increase of Jupiter's eccentricity. Also the secu- 
lar resonance g\ is stronger after the 2:1 MMR crossing for 
the higher eccentricity of Jupiter and also because of the 
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